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P q in small 

bin from its Poisson counterpart p q is studied within the 
Ginzburg-Landau description for second-order quark-hadron 
phase transition. Dynamical factor d q = P q /p q for the distri- 



bution and ratio D q 
behaviors between D, 



= dq/di are defined, and novel scaling 
are found which can be used to detect 



the formation of quark-gluon plasma. The study of d q and 
D q is also very interesting for other multiparticle production 
processes without phase transition. 

PACS number(s): 05.70.Fh, 05.40. +j, 12.38. Mh 



Multiplicity distribution is one of the most important 
and most easily accessible experimental quantities in high 
energy leptonic and hadronic collisions. From the well- 
known KNO scaling and its violation [1, 2] to the novel 
scaling form [3] investigated very recently the distribu- 
tion shows a lot about the dynamical features for the pro- 
cesses. Local multiplicity distributions have been stud- 
ied for many years in terms of a variety of phase space 
variables [4], and substantial progress has been made re- 
cently in deriving analytical QCD predictions for those 
observables [5]. Based on assuming the validity of the 
local parton hadron duality hypothesis, those analytical 
predictions for the parton level can be compared to ex- 
perimental data. A global and local study of multiplic- 
ity fluctuations [6] shows, however, that the theoretical 
predictions have significant deviation from experimental 
data. The significant deviation of theoretical predictions 
from experimental data indicates that we know only a 
little about multiparticle production processes since the 
hadronization process in soft QCD is far from being un- 
derstood. 

In this paper we try to investigate multiplicity distribu- 
tion in some small two-dimensional kinetic region (which 
can be rapidity and transverse momentum or azimuthal 
angle, for example) in high energy heavy-ion collision pro- 
cesses. In such collisions a new matter form, quark-gluon 
plasma (QGP), may be formed which subsequently cools 
and decays into the experimentally observed hadrons, 
thus the system undergoes a quark-hadron phase tran- 
sition. The hadrons produced in such processes may, in 
principle, carry some relic information about their par- 
ent state. Thus the investigation of the multiplicity dis- 
tribution may be interesting and useful for probing the 
formation of QGP. In this paper we are limited to dis- 
cussing multiplicity distribution under the assumption 



of a second-order phase transition, within the Ginzburg- 
Landau description for the phase transition. Within the 
same description for quark-hadron phase transition, the 
scaled factorial moments are studied by a lot of authors 
for second-order [7, 8] and first-order [9-12] phase transi- 
tions, and a universal scaling exponent v ~ 1.30 is given 
in [7,8,11,12]. 

It is useful to point out that the study of multiplicity 
fluctuations in photon production at the threshold of las- 
ing, which shows a similar type of phase transition [13], 
is already in its mature age, although the theory and 
experiment for a quark-hadron phase transition are still 
in their infancy. As explained in [7-12] the Ginzburg- 
Landau model, which has been used in describing su- 
perconducting transition and other macroscopic second- 
order phase transitions, can also be used to describe the 
multiplicity fluctuations in both second- and first-order 
phase transitions. In [12] multiplicity distributions are 
studied for both first- and second-order phase transitions. 
The authors showed that for second-order phase transi- 
tion the probability P q of finding q hadrons in the small 
bin under investigation decreases monotonically with the 
increase of q regardless of value of bin width and that for 
first-order phase transition P q is a decreasing function of 
q for small bin width whereas the shape of the distribu- 
tion changes with the increase of bin width. Thus the 
shape of the distribution was claimed a tool for telling 
the order of the phase transition. 

In this paper we first show that the criterion in [12] 
based on the shape of multiplicity distribution for the 
order of the phase transition is equivocal. This is eas- 
ily seen once one considers the trivial case without dy- 
namical fluctuations. For such a case, the multiplicity 
distribution p q is a Poisson one 
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■exp( 



(1) 



with s the mean multiplicity. From this distribution one 
has 



Pq+l 
Pq 



g + i 



(2) 



If s < 2.0 p q is a monotonically decreasing function of q 
whereas p q changes its shape for s > 2.0. Using the same 
parameters as in [12] s is calculated and listed in Table 1. 



Thus one can see that the shapes of multiplicity distri- 
butions given in [12] are similar to those of Poisson ones. 
In real experiments, one can always choose bin width 
to ensure the mean multiplicity larger than 2.0, then one 
cannot tell whether the distribution is shaped due to sta- 
tistical fluctuations or due to the dynamics in the phase 
transitions. So one cannot give the order of quark-hadron 
phase transition just from the general shape of the distri- 
butions, and detailed information is needed. This result 
is not surprising, because the non-dominant dynamical 
fluctuations can only modify the shape of statistical dis- 
tribution to some extent but cannot change its general 
behavior drastically. 

Nevertheless, it should be pointed out that the study of 
multiplicity distribution is still very interesting and use- 
ful for processes with the onset of quark-hadron phase 
transitions. In Ginzburg-Landau theory, the multiplicity 
distribution turns out to be a Poisson one if the field is 
purely coherent. Conversely, the distribution turns into 
a negative binomial one if the field is completely chaotic. 
In reality, one can assume multiplicity production aris- 
ing from a mixture of chaotic and coherent fields, so the 
multiplicity distribution in real processes is not a Poisson 
one nor a negative binomial one, and the deviation of the 
distribution from a Poisson one is due to dynamical fluc- 
tuations. The real quantity concerned is the deviation of 
experimental P q from its theoretical Poisson counterpart 
p q . Thus studying the deviation may reveal features of 
dynamical mechanism involved. Let the probability of 
having q hadrons in a certain bin is P qi the deviation of 
P q from its Poisson counterpart p q can be measured by a 
ratio d q = P q /p q . For the definition of d q to make sense, 
it is necessary to let the mean multiplicity s for P q and 
p q be the same. Dynamical fluctuations are shown to be 
existent if the ratio is far from 1.0, cither much larger 
or much smaller, for some q. The ratio d q can be called 
dynamical factor, since it is 1.0 unless there are dynam- 
ical fluctuations in the process. In the Ginzburg-Landau 
description for second-order phase transition P q is given 
by [8] 



P q (S) = Z- 1 / VfaqiS 2 I | 2 )e 



-F\, 



(3) 



where Z = J V(f>e^ F ^ the partition function, p q (s) the 
Poisson distribution with average s, and F[<f>] the free 
energy functional 

F[4>] = S 2 [a | | 2 +b | | 4 ] . 



It is instructive to note that a free energy functional 
with O(N) QCD order parameter was studied in [14]. 
The free energy functional used here is different from 
that in [14] because of the consideration that we are now 
only interested in the final state charged hadrons (most 



of them are n ) so that a two-component order param- 
eter is enough (which is written as a complex number) 
for our purpose. One can see that the functional used 
here can be derived from that in [14] by integrating out 
all other components and neglecting higher order powers 
of left components in the exponential. One more simpli- 
fication used in present functional is that the derivative 
term is neglected since former studies (see the last two 
papers in [7] for details) find out that the term has little 
contribution to the universal exponent which is a mea- 
sure of the fluctuations. Because of this simplification, 
the non-Gaussian functional integral can be treated as a 
normal integral and can be evaluated directly 

With the free energy functional above the system is in 
the plasma state for a > (the order parameter | <j>o | 2 
corresponding to the minimum of F[(j)\ is zero) and in 
hadron phase for a < (the order parameter | 0o | 2 > 0). 
In real experiments the temperature at which hadrons are 
emitted from the source is unknown and may be different 
from event to event. So we treat a as a free parameter and 
discuss only for a < in the following since in the quark 
phase with a > only a few hadrons can be produced 
through fluctuations. From the distribution of Eq. (3) 
one gets the mean multiplicity for a < 



_ Ji(| a | x) 
J (| a | x) 



(4) 



with J„(a) = j °° dyy"exp(-y 2 + ay), 



= 5/ ' \fb rep- 
resenting the bin width <5, a ex T — T c representing the 
temperature when the phase transition takes place. For 
small phase space bin the mean multiplicity in the bin is 
proportional to x thus can be very small. Under such cir- 
cumstance the distribution must be concentrated around 
Pq, and both P q and p q for q > 1 must be very small, 
so a direct comparison between them could induce large 
uncertainty. This demands that the bin width in real ex- 
perimental analysis should be large enough to ensure the 
mean multiplicity in the bin not too small (larger than 
0.5, say). Of course, smaller bins can be used if a precise 
determination of both P q and p q can be obtained from 
high statistical experimental data. For zero bin width the 
relevant results are rather sensitive to the cascading pro- 
duction of particles through resonances, so an extremely 
small bin width should be avoided. 

Because of the normalization of both P q and p q , the 
dynamical factor d q must be larger than 1.0 for some q 
and less than 1.0 for some other q if there exist dynamical 
fluctuations. One can easily derive 



d q (x) 



J q (x(\ a | -1)) J q (x | a |) 
J (x | a |) Jf(x | a |) 



exp(s) 



(5) 



The dependence of d q on q for different x and | a \ 
is shown in Fig. 1 by choosing | a |=1.0 and 2.0, 
— \nx = —1.0, 0.0, 1.0, 2.0, 3.0, respectively. From this 



figure, one can see that the general shapes of d q are sim- 
ilar for different choices of | a | but depend strongly on 
the bin width x. In detail, for large x (small —In a; or 
high mean multiplicity) d\ is quite large while d q> \ are 
smaller than 1.0. For small x (large —In a; or low mean 
multiplicity), however, d\ is smaller than but close to 1.0 
while d q> \ are larger than 1.0, indicating that two or 
more particles are more likely to be in the same small 
bin than for the pure statistical case. This phenomenon 
may be associated with the cluster effect or mini-jets in 
quark-hadron phase transition. For small x the values 
for d q are independent of | a |. From Fig. 1 one can also 
see that the dependence of d q on x is quite complicated. 
For some large q d q is monotonically increasing with the 
decrease of x, but for small q, d q first decreases and then 
increases with the decrease of x. Complicated behaviors 
can be seen for p q (x), considering the fact that s is an 
increasing function of x while p q (s) changes its behav- 
ior at s = q. But, the ratio p q {x)/p\{x) is an increasing 
function of x for q > 1. Moreover there exists a scaling 
law between p q 



p q (s) _ 21- 1 fp 2 (s) 



Pi( s ) 



^W 9-1 



PlO) 
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Thus it may be more interesting to study the dependence 
on x of 



D oS $ = 



PJPi 



dl Pq/Pl 



(7) 



instead of d q , and one may expect some scaling behav- 
ior of D q when the resolution is changed. Now we turn 
to study D q for second-order quark-hadron phase transi- 
tion. If there is no dynamical reason, P q = p q , one can 
see that D q for all q can have only one value, 1.0, no 
matter how large or small the bin width is. So from the 
range of values D q takes, one can evaluate the strength 
of dynamical fluctuations. D q can be expressed in terms 
of J n (ct) as 



lnD q = (g-l)ln 



Jo (| a\x) J q (x(\ a | -1)) 



hi 



M\a\x) Ji(x{\ a | -1)) ' 



(8) 



Besides x, there is in last expression another parame- 
ter | a | which is a measure of how far from the criti- 
cal temperature the hadronization process occurs and is 
unknown in current experiments. First let us fix | a \ 
to be 1.0. One can immediately see that, for any x, 
D q oc D^~ , a power law satisfied by Poisson, binomial 
and many other distributions in statistics. For any value 
of | a | D q increases monotonically with the decrease of 
x. This shows that the dynamical influence can be ob- 
served more easily in high resolution analysis. This can 
be understood physically, since different dynamical fluc- 
tuations may offset each other and become less obviously 
observable in large bin analysis. The behaviors of h\D q 



as functions of resolution — In x are shown in Fig. 2 for 
| a |=1.0 and 2.0 for g=2,3,4,5,6. For small — h\x val- 
ues of D q depend strongly on parameter | a |, but they 
approach parameter | a | independent values for large 
— In x. The similarity in the shapes of In D q as functions 
of — In x suggests a power law for other | a | between 
D q and D 2 similar to the case we showed for | o |= 1.0. 
In D q are rcshown in Fig. 3 as functions of In D 2 with 
the same data as in Fig. 2. For both | a |=1.0 and 2.0 
perfect linear dependences of hiD q on \1\D2 can be seen. 
For other values of | a | the similar linear dependence is 
checked to be true. Thus one has 



\nD q = A q 



B q \nD 2 



(9) 



with A q and B q depending on | a |. The fitted results 
of A q and B q from curves in Fig. 3 are shown in Fig. 
4 as functions of ln(g — 1) for | a |=1.0 and 2.0. It is 
obvious that both In A q and In B q have linear dependence 
on ln(o — 1) for fixed \ a \. Especially, for the purpose of 
studying power law, we investigate B q and find that 



B„ 



(?-l) 7 , 



(10) 



with 7 depending on | a |. For visualization, the linear 
fitting curves for In B q vs ln(g — 1) are shown also in Fig. 
4 for I a | = 1.0 and 2.0. The slopes for \nA q is about 
twice those for hiB qi and they increase with increasing 
I a |. When | a | is zero, corresponding to the case in 
which hadrons are produced exactly at the critical point, 
numerical results show that 7 is 0.819. With the increase 
of I a I, 7 increases quickly. For sufficiently large | a |, 
when the difference between | a \ — 1 and | a | can be 
neglected, corresponding to the case in which hadrons are 
produced at temperature much below the critical point, 
one finds that D q (x) = F q (\ a \ x), with F q the scaled 
factorial moment which is given in [8] for second-order 
phase transition as 

F q (x)^J q {x)Jl 1 (x)J-\x). 

Similar relation between D q and F q is also true in the 
small x limit. In these limiting cases, the scaling of D q 
is equivalent to that of the scaled factorial moments F q , 
and one can get the exponent 7=1.3424 for large — ln.x 
[15] or large \ a \. The dependence of 7 on | a \ is shown 
in Fig. 5. In real experiments, | a \ is not known for an 
event and may be increasing in the hadronization process. 
Thus some average over | a \ should be made. The smaller 
I a I, the less the number of produced particles. Thus the 
main contribution to P q comes from events with large 
I a I or with high multiplicities. For those events, one 
should get 7 near 1.30, close to the universal exponent 
v given in former studies of F q for second-order phase 
transition. For events with low multiplicity, one can get 
7 > 0.819. So the theoretical range for the exponent 7 
is (0.819, 1.3424), corresponding to temperature range 
from T = T C toT<T c . 



In conclusion, two new quantities d q and D q are in- 
troduced to describe dynamical fluctuations in quark- 
hadron phase transitions. In the Ginzburg-Landau de- 
scription for second-order quark-hadron phase transi- 
tion, d q and D q are investigated analytically, and it is 
found that D q obeys a power law, D q ex D 2 q , with 
B q = (q — l) 7 . In experimental analysis, both d q and D q 
can be obtained quite easily. To get P q one only needs to 
count up the number of events with exact q hadrons in 
the bin. p q is of Poisson type and can be calculated from 
the experimental s. Simple algebras give d q and D q . The 
existence of dynamical fluctuations can be confirmed if 
d q and D q can take values very different from 1.0. The 
scaling between D q and D2 is a possible signal for the 
formation of QGP, because up to now no other dynami- 
cal reason is known to induce such a scaling. The value 
of exponent 7 can be used to measure the deviation of 
the temperature, at which the hadronization occurs, from 
the critical point. The study of D q should be carried out 
in real experimental analysis in the future to see whether 
QGP has been formed in current high energy heavy-ion 
collisions. As a tool to study dynamical fluctuations d q 
and D q introduced in this paper may also be interesting 
in experimental analysis of leptonic and hadronic inter- 
actions without quark-hadron phase transitions. Study 
of dq and D q in first-order quark-hadron phase transition 
is in preparation. 
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TABLE 
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1 


2 


4 


6 


B=+l 


0.226 


0.429 


0.804 


1.155 


B = -l 


0.342 


0.797 


2.172 


4.582 



TABLE I. Mean multiplicities s for second-order (B—+1) 
and first-order (B — —1) phase transitions for different bin 
widths, x is a parameter (different from the quantity used in 
this paper) associated with bin width 8, parameter s in [12]. 
The mean multiplicities are calculated using Eq. (13) in [12]. 
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Figure Captions 

Fig. 1 Dependence of dynamical factor d q on q for 
| a |=1.0 and 2.0, for -In x = -1.0,0.0,1.0,2.0,3.0. 

Fig. 2 Dependences of D q on bin width — lnx for | a | = 
1.0 and 2.0 for q=2, 3, 4, 5, 6. 

Fig. 3 Scaling behaviors between D q and D^ for | 
a | = 1.0 and 2.0. The data are the same as in Fig. 
2. From lower to upper are curves for <?=3, 4, 5, 6, 
respectively. 

Fig. 4 Coefficients for the scaling between D q and Z?2, 
In D q = A q + B q In D2, as functions of \n(q — 1) for 
I a |=1.0 and 2.0. Linear fitting curves are shown 
for B q = {q - 1)T. 

Fig. 5 Dependence of exponent 7 on | a |. For large 
I a J, 7 is about 1.34. 
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